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Abstract 

An established procedure for the estimation of ice-water interfacial energy based on evaluation of the slope 
of the experimental ice nucleation rate data versus scaled temperature is critically analyzed in this work. An 
inconsistent estimate of the ice-water interfacial energy is found in the work of Murray et al. [Phys. Cliem. 
Clrem. Phys., 2010, 12, 10380-10387]. The source of the inconsistency is identified in an inappropriate 
regression method used for experimental ice nucleation data fitting, a correct estimate of the ice-water 
interfacial energy is presented, and limits of the slope method are discussed. 


1. Introduction 

The classical nucleation theory (CNT) provides 
a theoretical link between the ice nucleation rates 
and the thermophysical properties of the ice-water 
system, i.e. the density of the solid phase, the diffu¬ 
sion coefficient of the liquid phase, the equilibrium 
pressure, and the interfacial energy [6]. Therefore, 
in cases where experimental measurements of the 
nucleation rate are available, CNT presents an in¬ 
direct, theoretical approach for the estimation of 
the interfacial energy provided that the other above 
mentioned thermophysical properties are known. 

For ice and water, the interfacial energy mea¬ 
surements reported in literature are restricted to 
the triple point temperature [3], The estimation 
of the interfacial energy for lower temperatures in 
the supercooled region, was possible only due to 
ice nucleation rate measurements HU or molecular 
simulations of ice nuclei formation Pla¬ 
in this work, a method for the estimation of ice- 
water interfacial energy from the experimental nu¬ 
cleation rates mm will be investigated that is 
based on the evaluation of the slope of the loga¬ 
rithm of the measured nucleation rate In J vs. the 
function T -3 (ln(S 1 )) -2 of the experimental temper¬ 
ature T, where S is the ratio of water and ice satu¬ 
ration pressures. It will be referred to as the slope 


method in the following text. Several authors es¬ 
timated the ice-water interfacial energies based on 
the slope method nn Ha gb iu covering the temper¬ 
ature range 200 240 K. I will study the validity 

of the ice-water interfacial energy estimates by the 
slope method in this paper. Particularly, I intend to 
show that the ice-water interfacial energy estimate 
of Murray et al. HU suffers from internal inconsis¬ 
tency and needs to be reconsidered. The nature of 
this inconsistency will be studied in detail to point 
out the weak spots of the slope method. 

The paper is arranged in the following way. First, 
the equations of the slope method for interfacial en¬ 
ergy estimation are derived in section [2] In section 
[3j a proof is presented showing the inconsistency of 
the interfacial energy estimate by Murray et al. HU- 
In section [4| the source of the inconsistency is iden¬ 
tified, a correct estimate of the ice-water interfacial 
energy is presented, and the limits of applicability 
of the slope method are discussed. 


2. Slope method 

The slope method mm for the estimation of 
the ice-water interfacial energy is based on the CNT 
nucleation rate equation [7] 
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J = Jo exp 
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~kT 
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relating the nucleation rate J [m 3 s x ] to the nu- 


cleation work of the critical cluster W* [J]. The pre¬ 
factor J 0 [m _ 3 s _1 ] reflects the kinetics of the cluster 
growth and T [K] is absolute temperature. By us¬ 
ing the CNT expression for the critical nucleation 
work W* = 3 ( 1 fcj] n g) 2 ; Eq. (1) can be rearranged 
to 

167tj 3 v 2 


In J = In Jn — 


3k 3 T 3 (lnS) 2 


( 2 ) 


where 7 [J/nr 2 ] is the ice-water interfacial energy, v 
[nr 3 ] is the molecular volume of the solid phase, and 
the ratio of saturation pressures S = p{ q /p e s q rep¬ 
resents the supersaturation of the liquid. Here, the 
vapor-liquid equilibrium pressure is denoted p/ 9 , 
and the vapor-solid equilibrium pressure is denoted 
p eq . The vapor-solid equilibrium pressure of cubic 
ice m is used in this work to retain consistency 
with the work of Murray et al. PH- However, the 
structure of the clusters in ice nucleation is still un¬ 
der debate ITU] . 

To transform Eq. © into the relations of the 
slope method, the following assumptions are made. 
The interfacial energy is taken constant, i.e. in¬ 
dependent of temperature. The density of the solid 
phase is taken constant. And the pre-factor is taken 
constant. Since the size of the temperature inter¬ 
val of the experimental nucleation data analyzed by 
the slope method is typically a few Kelvins, the as¬ 
sumptions of the constancy of the thermophysical 
properties over this narrow temperature interval are 
plausible, and Eq. © can be written in the form 

In J = n + mt s (3) 


where the parameters 


n = In Jq 


(4) 


and 


to = — 


167T7 3 u 2 

3 k 3 


(5) 


are both due to the above assumptions constant, 
and 


1 

T 3 (ln S) 2 


( 6 ) 


is a scaled temperature. As a result, the slope 
method is simply a fit of the experimental ice nu¬ 
cleation rate data to the linear function © in t s , 
giving the parameters n and to. The ice-water in¬ 
terfacial energy 7 corresponding to the analyzed ice 
nucleation data is recovered from the slope to ac¬ 
cording to Eq. ([ 5 ]), and the pre-factor Jq is given by 
the absolute value n according to Eq. ©• 



T [K] 


Figure 1: Experimental ice nucleation rates reported by 
Murray et al. El and Stan et al. US- The full lines show 
fits of the experimental data, and the dashed lines delimit 
the standard deviation bands. The fit of Murray et al. El 
is used as presented in their Fig. 4. The data of Stan et 
al. ES as presented in their Fig. 8 are shown with the re¬ 
ported ±0.4 °C standard deviation. 


3. Inconsistency of interfacial energy esti¬ 
mates 

The validity of the slope method estimates of 
the ice-water interfacial energy based on the exper¬ 
imental nucleation data as published by Murray et 
al. m will be investigated in this section. 

I have identified inconsistent estimates of inter¬ 
facial energy for two specific sets of ice nucleation 
data investigated by Murray et al. Cl, i.e. their 
own experimental data and the data measured by 
Stan et al. Il8j , respectively. The two experimental 
data sets are plotted in Fig. [T] By using the slope 
method, Murray et al. m estimated ice-water in¬ 
terfacial energy 7 m = 20.8 ±1.2 mJ/m 2 from their 
own experimental data, and ice-water interfacial en¬ 
ergy 7 s = 23.7 ± 1.1 mJ/m 2 from the experimental 
data of Stan et al. [15] , respectively. It is the gap 
between the two intervals, 7 m and 7 s, that renders 
the two interfacial energy estimates inconsistent. In 
other words, as I will show below in detail, since the 
initial two experimental datasets used as input for 
the evaluation of the respective interfacial energies 
do overlap to a large extent, which is seen in Fig. 
[l] it is impossible to calculate two interfacial en¬ 
ergy estimates by the slope method that form two 
disjunct intervals, as is the case with and 75 . 

To prove the particular inconsistency in interfa¬ 
cial energy estimates 7 m and 7 g suggested above, 
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I will frame general requirements for a physically- 
relevant deduction of the ice-water interfacial en¬ 
ergy from the nucleation rate according to the CNT. 
Mathematically, the deductions will utilize elemen¬ 
tary properties of continuous, strictly monotonic 
functions in a reductio ad absurdum type of log¬ 
ical proof. The slope method analysis will present 
a limiting case of this general scenario. 

Let us denote g a continuous, strictly monotonic 
function of temperature g: T —» 7, which stands for 
the functional dependence of the interfacial energy 
on temperature. Similarly, let j be a continuous, 
strictly monotonic function of temperature j: T —7 
J, which stands for the functional dependence of 
the nucleation rate on temperature, e.g. in the form 
of Eq. ([I]). Let the temperature interval [Ti, T 2 \ be 
a subdomain of the two functions g and j that corre¬ 
sponds to the temperature range of the analyzed ice 
nucleation experiment. Both functions g and j are 
bijective by their above definitions. In other words, 
they present a one-to-one correspondence between 
the elements of subdomain [Ti, T 2 ] and the elements 
of functional images g([Ti,T 2 ]) and j([Ti,T 2 ]), re¬ 
spectively. Therefore, inverse function j -1 : J —> 
T exists and it is bijective as well. Furthermore, 
composite function (jj = g(j~ 1 )'. J —> 7 is bijective. 
Here, gj is the formal representation of a general 
theoretical method to deduce the interfacial energy 
from the nucleation rate in temperature range pi, 
Ti]- 

Now, let = [J n , J 12 ] and J 2 = [J 2 i, J22] be 
two intervals of nucleation rates that possess a non¬ 
empty intersection Jo = JiO J 2 7^ 0 . Let 71 = [711, 
712] and 72 = [721, 722] be gq’s functional images 
of intervals J 1 and J 2 , i.e. 71 = gj(Ji) and 72 
= 5 j(J 2 ), respectively. Intervals 71 and 72 repre¬ 
sent the ranges of interfacial energies deduced from 
the two nucleation rate data sets J\ and J 2 ac¬ 
cording to a general method gj. Then VJ £ Jo it 
holds that gj (J) £ 71 and gj(J) £ 72- In other 
words, for any value of the nucleation rate J, which 
belongs to the intersection of intervals Ji and J 2 , 
the deduced value of interfacial energy according to 
method g 3 belongs to interval 71 and also to interval 
72. Therefore, since element 7 = gj (J) belongs to 
both intervals 71 and 72, intersection 70 = 71 H 72 
has at least one element 7 and, therefore, 70 must 
be non-empty, i.e. 70 7^ 0. 

For the particular case of the ice nucleation ex¬ 
perimental data of Murray et al. EH and Stan et 
al. [ 18 ] the intersection of the measured nucleation 
rate data sets is not empty, as can be seen in Fig. |T] 


the measured values of the nucleation rate in the 
two datasets clearly overlap to a large extent. Ac¬ 
cording to the above reasoning the range of interfa¬ 
cial energies deduced from the Murray et al. m nu¬ 
cleation data must have an non-empty intersection 
with the range of interfacial energies deduced from 
the Stan et al. [lS] nucleation data provided that 
the method used to deduce the interfacial energy 
assures that the temperature dependency of the in¬ 
terfacial energy and the nucleation rate is strictly 
monotonic. The strict monotonicity of both the 
nucleation rate and the interfacial energy is a natu¬ 
ral requirement conforming to the physical reality; 
the nucleation rate decreases with increasing tem¬ 
perature and the interfacial energy increases with 
increasing temperature [ 8 ]. 

However, the slope method assumes, rather 
unphysically, that the interfacial energy is con¬ 
stant, and not strictly monotonic as supposed in 
the general case above. Under such assumption, 
the considerations presented above result into an 
even stronger requirement on the deduced interfa¬ 
cial energies from two intersecting nucleation rate 
datasets. Any constant function can be viewed as 
a limit of a sequence of strictly monotonic func¬ 
tions, e.g. a sequence of linear functions f n = 
C + a n (T — To) with slope a n decreasing as 1/n 
for n —> 00 . For every element of the sequence 
the above reasoning applies and the deduced in¬ 
terfacial energy ranges must possess a non-empty 
intersection. In the limit of the constant function 
the intersection shrinks to a single value of interfa¬ 
cial energy, common to both datasets. As a result, 
the slope method must produce interfacial energies 
equal to each other within their standard deviation 
ranges from two overlapping datasets of nucleation 
rates. However, we already stated that the esti¬ 
mates of the interfacial energy of the slope method, 
7 m and 7 s, as published by Murray et al. m are 
not equal within their standard deviations, which 
is in contradiction with the outcome of the above 
reasoning. Therefore, the slope method procedure 
as preformed by Murray et al. El] must contain 
a hidden inconsistency, an erroneous step that in¬ 
troduces additional uncertainty. Its nature will be 
discussed in the following section. 

4. Discussion 

The fit of the experimental ice nucleation-rate 
data to the linear function © forms the basis of 
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Figure 2: Comparison of the ordinary least squares regres¬ 
sion method (dashed line) and the orthogonal distance re¬ 
gression method (full line) for the estimation of the slope of 
the ice nucleation data of Murray et al. Bli¬ 
the slope method. The fit results in two parame¬ 
ters n and to that will be discussed in this section. 
First, the evaluation of the slope m resulting in the 
estimation of ice-water interfacial energy © will be 
investigated with the goal of resolving the inconsis¬ 
tency found in Sec. [3j Second, the evaluation of the 
nucleation rate pre-factor from the term n accord¬ 
ing to Eq. Q will be analyzed. 

4-1. Interfacial energy 

Although Murray et al. HU do not specify the 
particular regression method that they used to per¬ 
form the linear fit to their data yielding their re¬ 
ported slope moLS = — (6.02 ± 0.36)xl0 7 K 3 , 
their result can be reproduced with the widely used 
ordinary least squares (OLS) method. By using 
the OLS method implemented in Python (pack¬ 
age scipy. stats . linregress) to fit the Murray 
et al. E! experimental data to the function ([3]) we 
can check that the above mentioned value of tools 
is recovered. Therefore, let us assume that the OLS 
method was indeed used by Murray et al. El in 
their analysis. 

The OLS method was derived with the assump¬ 
tion that the observation of the independent vari¬ 
able, i.e. the temperature in this case, is error-free. 
From this point of view the applicability of the OLS 
method is highly questionable for nucleation-rate 
experimental data, because the error in the temper¬ 
ature measurement of an ice nucleation experiment 


is typically larger than ±0.4 °C, and it is identified 
as the main source of uncertainty in the nucleation 
measurement [18] . 

The regression method derived with the assump¬ 
tion of non-zero observational errors in the indepen¬ 
dent variable is the orthogonal distance regression 
method (ODR), also known as errors-in-variables 
modeling, or total least squares [ 5 j. The Python 
ODR implementation (scipy.odr.odrpack) ap¬ 
plied to Murray et al. El ice nucleation data results 
in the slope mo dr = — (8.04 ± 0.48) xlO 7 K 3 . 

The difference in the slope estimates between the 
OLS and ODR methods for the Murray et al. El 
ice nucleation data is shown in Fig. [2] The ice- 
water interfacial energy corresponding to the moDR 
is 7 odr = 22.9 ± 0.5 mJ/m 2 according to the slope 
method, Eq. ([5]). And after including the uncer¬ 
tainty in the cubic ice sublimation pressure |17| (± 
0.8 mJ/m 2 ), the final estimate of ice-water interfa¬ 
cial energy using the ODR method is 7 odr = 22.9 
± 1.3 mJ/m 2 . This new ODR estimate is roughly 
10 % higher that the OLS estimate 7 ols = 20.8 ± 
1.2 mJ/m 2 calculated from tools- The new value 
70 dr is sufficiently close to the value 7 s Murray et 
al. El estimated from the nucleation data of Stan 
et al. M > and therefore 7 o dr does not suffer from 
the contradiction described in Sec. 0 

Note that the Stan data HE! are fitted correctly 
even with the OLS method that fails for fitting the 
Murray data El- We can check that the OLS 
method gives the slope m = —8.952xl0 7 K 3 and 
the ODR method results in m = —8.967xl0 7 K 3 , 
which both correspond to the interfacial energy 75 
= 23.7 mJ/m 2 reported by Murray et al. |TT] for the 
nucleation data of Stan et al. [T5 . The Stan data, 
as presented in Fig. [I] (and in their Fig. 8 [TB]h on 
the contrary to the raw data by Murray et al. El, 
are already averaged from a set of more than 37 
thousand freezing experiments, effectively removing 
the large uncertainty in experimental temperature 
measurement. The Stan data are therefore suitable 
for fitting even by the OLS method. 

4-2. Nucleation rate pre-factor 

The slope method concurrently fits two parame¬ 
ters, n and to, from a given data set of experimental 
nucleation rates. The two parameters are therefore 
coupled and we can expect that the larger an er¬ 
ror in the estimation of to we make, the larger the 
corresponding error in n we get. It was shown in 
the preceding section that the interfacial energy es¬ 
timate 7 m suffered from an error due to the inap- 
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Figure 3: Pre-factor Jo estimated by Murray et al. ED with 
the slope method (as given in their Table 2) from the ice 
nucleation data available in literature (red circles). Horizon¬ 
tal errorbars show the temperature range of the input ex¬ 
perimental nucleation-rate datasets. The vertical errorbars 
denote the uncertainty of Jo fit. The green full line shows 
the pre-factor according to the volume-diffusion controlled 
growth model [7]. The blue square shows the OLS-based 
fit of the pre-factor, and the the green diamond shows the 
ODR-based fit of the pre-factor from the experimental data 
of Murray et al. im¬ 
propriate statistical processing of the experimental 
data. Therefore, an error was introduced in the cor¬ 
responding estimate of nucleation rate pre-factor J 0 
as well. 

The difference between the pre-factor estimates 
using the OLS and the ODR regression methods is 
shown in Fig.[3]for the experimental data of Murray 
et al. HO- The fitted parameters are hols = 75 ± 
5, and no dr = 94.5 ± 5, respectively. By evaluat¬ 
ing the pre-factor according to Eq. Q , we find that 
the OLS-based estimate of Jo is roughly 9 orders of 
magnitude lower than the ODR-based result. It is 
the ODR-based estimate of Jo that is within its 
uncertainty equal to an independent model of Jo, 
i.e. the volume-diffusion based description of the 
cluster growth [7] (p. 141), as shown in Fig. |3j 

4.3. Final remarks 

Not only the already discussed regression issue is 
the source of uncertainty in the slope method. The 
assumptions of the method itself, as summarized in 
Sec. [2j play a role as well. The constancy of the 
physical properties assumed by the slope method is 
plausible for small temperature ranges of the ana¬ 
lyzed experimental nucleation rates. As the tem¬ 


perature range of the experiment increases, the in¬ 
terfacial energy, ice density, and the pre-factor are 
getting less accurately approximated by a constant. 
It is therefore desirable for the slope method to an¬ 
alyze temperature ranges as small as possible. On 
the other hand, to evaluate the slope of the ice nu¬ 
cleation data from the cloud of scattered experi¬ 
mental nucleation rates one needs as large temper¬ 
ature ranges as possible, because the uncertainty 
of the fit of the slope gets larger as the temper¬ 
ature range of the analyzed data decreases. Ob¬ 
viously, both requirements, i.e. having the tem¬ 
perature interval small for the assumptions of the 
slope method to be valid and having the temper¬ 
ature interval large for a precise evaluation of the 
slope, cannot be simultaneously satisfied. And un¬ 
certainty due to this fact will be inevitably present 
in the results of the slope method. 

Fig-i shows the Jo estimates from all the ice 
nucleation-rate datasets in the temperature range 
235 - 240 K available in literature as calculated 
by Murray et al. HJ. A huge scatter in the pre¬ 
factor estimates spanning roughly 25 orders of mag¬ 
nitude can be observed. However, the pre-factor 
is only slightly temperature dependent in the case 
of nucleation in liquids [7 (p. 199). According to 
the above-mentioned diffusion-based growth model 
the pre-factor depends linearly on the self-diffusion 
coefficient of supercooled water, which was mea¬ 
sured down to 238 I< with accuracy less than 1% 
m- Therefore, the volume-diffusion based growth 
model gives us an estimate of the pre-factor Jq with 
negligible uncertainty compared to the uncertainty 
of 25 orders of magnitude reported in Murray et 
al. im- Obviously, the scatter of Jo estimates as 
derived by Murray et al. m and shown in Fig. [3] is 
non-physical, and reflects the errors in implement¬ 
ing the slope method. 

Finally, the slope method was also used by Manka 
et al. [9] to estimate the ice-water energy of 15.6 
mJ/m 2 from their ice nucleation data measured in 
the temperature range 202 - 215 K. The authors 
do not report any uncertainty of their interfacial 
energy estimate; they just state ”we analyzed all of 
our data using the same formalism” as Murray et 
al. m • Therefore, a similar regression issue as in 
the above-discussed case of Murray et al. Cl inter¬ 
facial energy estimate is likely to arise. Also, the 
relatively large temperature interval of the Manka 
data collides with the assumptions of constancy of 
the thermophysical properties of the slope method. 
But most importantly, the CNT formulation used 
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in the slope method as presented in Sec. [2] is not 
valid for ice nucleation data measured at such low 
temperatures due to the omitted pressure effect, as 
shown recently [13] . The pressure-related terms in 
the CNT formulation present a considerable con¬ 
tribution to the nucleation work and, therefore, 
account for a shift in the CNT-predicted ice nu¬ 
cleation rate at low temperatures. The errors in 
CNT formulation are then projected to errors in the 
ice-water interfacial energy estimates of the slope 
method. 

5. Conclusion 

A method for the estimation of the ice-water in¬ 
terfacial energy based on the classical nucleation 
theory and on the evaluation of the slope of the ex¬ 
perimental ice nucleation rates vs. a scaled temper¬ 
ature was thoroughly analyzed in this work. A logi¬ 
cal contradiction in the interfacial energy estimates 
was identified in the work of Murray et al. [TT| that 
is related to the linear regression algorithm utilized 
for the evaluation of the slope of experimental ice 
nucleation rate data. The contradiction was re¬ 
moved by using the orthogonal distance regression 
method for a proper evaluation of the slope instead 
of the ordinary least squares method. The corrected 
estimate of the interfacial energy is by 10 % (2.1 
mJ/rn 2 ) higher than the original value reported by 
Murray et al. m. 

In the light of the findings of this work, esti¬ 
mation methods mm utilizing the absolute val¬ 
ues of the experimental nucleation rates instead of 
the temperature derivative, and using a theoretical 
model of the nucleation rate pre-factor instead of 
its fit from experimental nucleation data, present 
a safer way of deducing the interfacial energy from 
nucleation rate data avoiding uncertainties inher¬ 
ently contained in the slope method. 
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